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Abstract: We construct an effective three dimensional theory for the MSSM at high temperatures 
in the hmit of large-m^. We analyse the two- loop effective potential of the 3D theory for the case 
of a light right handed stop to determine the precise region in the mh-mi^ plane for which the 
] sphaleron constraint for preservation of the baryon asymmetry is satisfied. We also compare with 
^ ! results previously obtained using 3D and 4D calculations of the effective potential. A two-stage 
- - i phase transition still persists for a small range of values of nni^. The allowed region requires a value 
of mi^ < mt and m?, < 100 (110) GeV for tuq = 300 GeV (1 TeV). 



^On leave of absence from the Universidad Antonio Narino, Santa Fe de Bogota, COLOMBIA. 



1 Introduction 



The analysis of the electroweak phase transition in the MSSM has been the subject of intensive 
study in recent years. The main motivation is the generation (and preservation) of a possible baryon 
asymmetry at the electroweak scale |jl|, |^ (for reviews, see Many different contributions 

have given a clearer idea as to where in parameter space the condition of a strong enough first-order 
phase transition can occur in order for electroweak baryogenesis to be possible. Mainly two different 
analytic approaches have been used in the analysis of the phase transition for the MSSM. The first 
one involves the evaluation of the effective potential in the 4D theory. The one- loop analysis @-0 
led to the conclusion that low values of the ratio of the vacuum expectation values of the two Higgs 
doublets tan f3 = and large values of the pseudoscalar mass niA were favoured. More recently the 
effects of a very light right handed stop with a soft SUSY-breaking mass value of < were shown 
to strongly strengthen the phase transition by enhancing the cubic term in the effective potential 
[ p!0| , pJ]] . As first pointed out by Espinosa ||12[, two-loop corrections were also shown to be very 
important in order to increment the allowed parameter space for which electroweak baryogenesis 
can take place. In particular two-loop QCD corrections from stops can strongly affect the value of 
the scalar field at the phase transition |T^-[Q. A precise determination of the region in parameter 
space for which electroweak baryogenesis is viable was done by Carena et al. |T3], using a two- 
loop calculation in 4D, with a light right handed stop, and a heavy third generation left squark 
doublet, which is decoupled from the thermal bath. They conclude that the Higgs mass must be 
lower than 105 GeV and that the right stop mass must be in the range of 110-160 GeV if absolute 



stability of the physical vacuum is required. A very recent paper by Cline and Moore [0 use the 
two-loop effective potential in 4D, fully incorporating squark and Higgs boson mixing, to determine 



the allowed region in parameter space. Their results are in good agreement with those of ref. |15 



In the second approach for the analysis of the phase transition, the aim is to separate the 
perturbative and non-perturbative aspects of the theory. The purely perturbative component of 
the calculation is performed by constructing an effective 3D theory for the model under consideration 
r^-[^. The parameters in the 3D Lagrangian are obtained using dimensional reduction at high 



temperature by matching the static Green's functions in the two theories, to a given order in the 
perturbative expansion, by integrating out the non-zero Matsubara modes with masses of the order 
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of ttT, where T is the temperature. A further reduction can also be performed noting that some 
of the static modes in the theory have acquired thermal masses proportional to a gauge coupling 
multiplied by the temperature, ~ QwT, QsT. These so-called heavy particles can then be integrated 
out as well. The effective potential calculated in the 3D theory reproduces the results obtained with 
the 4D effective potential. References give more details concerning the construction 

of effective theories for both the Standard Model and the MSSM. This approach simplifies the 
theory as it is now purely bosonic, which facilitates lattice simulations. Numerical studies of the 
reduced theory will take into account the non-perturbative effects. We refer the reader to the above 
publications for further discussion regarding strengths and weaknesses of 4D and 3D calculations. 

For the generic case in which there is a single light scalar doublet field at the phase transition 
the constraint on strength of the transition is translated into an upper bound on the ratio of the 3D 



Higgs self-coupling and the square 3D gauge coupling [21|. This ratio has a weak dependence 
on the temperature for values close to the critical temperature of the phase transition for the 
allowed range of values of the Higgs mass [^. In addition, the 3D scalar and gauge couplings are 
renormalization group invariant. This implies that a 1-loop matching of the 3D coupling constants 
and masses to the physical parameters and the temperature suffices to determine the strength of the 
phase transition, using the constraint given by the non-perturbative analysis of the phase transition 
for a single light scalar field ETf . 



To clarify an essential point, we recall that the critical temperature for the transition from the 
high temperature minimum to the standard electroweak minimum is obtained from the condition 
that the value of the effective potential at these two minima are equal, ^(0, T) = V{(pc, T). Therefore 
we insist that the value of the critical temperature does depend on a precise determination of the 3D 
mass parameter, which is not renormalization group invariant. The scalar mass requires ultraviolet 
renormalization and a 2-loop calculation (in 4D) must be performed even for the case of a single 
light scalar at the phase transition. However, since the ratio of the vacuum expectation value of 
the scalar field to the temperature, ^ (or equivalently, ^), which determines the rate of sphaleron 
transitions in the broken phase has only a weak dependence, in this case, on the temperature, the 
two-loop calculation is not necessary. 

In the initial analysis for the MSSM using the 3D approach |^4[-|^, the light stop scenario 
could only be investigated for values of the right handed stop soft SUSY-breaking mass of mu > 50 
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GeV. For lower values of the stop mass the perturbative procedure of integrating out the "heavy" 
modes starts to break down as the relevant expansion parameter is of the form In addition, the 
numerical constraint from non-perturbative studies is no longer valid and new studies that include 
the effect of the SU(3) gauge fields must be performed. 

An effective Lagrangian for a light Higgs and a light stop must be constructed to analyse the 
phase transition for lower values of mu within the 3D approach [Q. A surprising result of the 
perturbative analysis in this scenario was that a possible two-stage phase transition could take 
place, in which the Universe would first undergo a transition to a colour-breaking vacuum and, at 
a lower temperature, another transition to the physical vacuum would occur. In fact, the work of 
Bodeker et al. was the first to point out that, for a small range of values of rn^^, the phase transition 
could occur in two stages for the MSSM. This analysis was based on a 2-loop calculation of the 
effective potential in 3D including the leading corrections in the dimensional-reduction procedure. 

It is of course of great interest to exactly identify the range of parameter space for which this two- 
stage phase transition can occur. In this case it is necessary to determine the critical temperatures 
very precisely for the two possible directions of the transition. The most relevant quantities that 
determine the critical temperatures are the 3D mass parameters for the Higgs doublet and the right 
handed stop. These masses depend logarithmically on two parameters, Ah^ or Ajj^, which can only 
be determined exactly from the 4D theory. In the initial reduction implemented by []T1|, the exact 
values of Afj-^^ and Ajj.^ were not determined. An estimate was used instead, based on the value of 
the corresponding parameter in the Standard Model. Here we employ a combination of 4D and 3D 
calculations of the effective potential to obtain the exact values of Ah-^ and A^/g. 

Undoubtably, many of the questions that arise from the limitations of perturbation theory 
will only be answered when the corresponding lattice calculations are done. However, the first 
non-perturbative results show that the 2-loop results for the strength of the phase transition are 
conservative in the bounds they impose on the allowed masses for the Higgs and the light stop 
These results were obtained using a simplified model in the reduction procedure. For more 



complicated initial Lagrangians the final effective 3D theory is characterized by the same couplings 
and masses. We would like to point out a few of the features that can be affected in the perturbative 
calculation of the effective theory and their consequent effect on baryogenesis. Additionally, in order 
to apply the non-perturbative results, a precise mapping of the 4D to 3D parameters is needed. 
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which is presented in this paper. The results obtained with the perturbative 2-loop effective potential 
presented here can then be compared with the lattice analysis. 

The procedure of constructing an effective 3D theory is based on the mass hierarchy which ap- 
pears at finite temperature. The validity of the results depends on an adequate expansion parameter 
and the suppression of the higher order terms in the 3D theory. The value of the masses of the 
particles which are integrated out will define the regime of validity of the approach. As mentioned 
above, previous analyses ||2^-[13] 3D theory integrating out the right handed stop claimed 



that for mu > 50 GeV the higher order terms were suppressed and consequently the effective the- 
ory was an adequate description. We compare the results obtained for the ratio of the vacuum 
expectation value of the Higgs field to the temperature using the two-loop effective potential with 
and without integrating out the right handed stop to determine more precisely the value of the 
right handed stop mass for which the effective theory is no longer valid. We find that for values of 
mu < 123 GeV the results obtained having integrated out the right handed stop are unreliable. 

The paper is organized as follows: in section 2.1 we present the dimensional reduction to the 
effective bosonic theory at one-loop. Section 2.2 presents a further one-loop reduction in the 3D 
theory, eliminating the heavy fields. In section 2.3 we give the expression for the 2-loop unresummed 
effective potential in 4D, which is necessary for evaluating Kh^ and Kuz- The contribution to the 3D 
effective potential from the "heavy" particles that were integrated out at the second stage is given 
in section 2.4. Section 2.5 discusses the relevant zero-temperature effects that must be included in 
our analysis. In section 3 we present our results for the critical temperatures and the strength of 
the phase transition. The allowed region for electroweak baryogenesis to occur is also given here. 
A comparison of the results obtained from the effective potential with and without integrating out 
the right handed stop is given in this section. Finally, in section 4, we conclude. The appendix 
contains the relevant formulae for the case in which the right handed stop is integrated out. 

2 Dimensional Reduction 

We will now perform dimensional reduction by matching, as has been previously done in refs. 



21],|Q-|26| for different models. Our initial 4D Lagrangian corresponds to the MSSM in the large- 



niA limit. The particles that contribute to the thermal bath are the Standard Model particles plus 
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third-generation squarks: tL,bL,tji,bji. We will only consider here the case of zero squark mixing. 
The results for non-zero squark mixing will be presented elsewhere^ |2^. There are two stages of 
reduction. The first one corresponds to the integration out of all non-zero Matsubara modes, that 
is with a thermal mass of the order of ~ ttT. We calculate all one-loop contributions to mass terms 
and coupling constants of the static fields to order g^, where g denotes a gauge or top Yukawa 
coupling. The second stage of reduction corresponds to the integration of heavy particles with 
masses of the order of gwT, gsT. 

2.1 First Stage 

The potential in the 3D effective theory after integration over non-zero Matsubara modes is of the 
form 



+ m%.^Q^Q + ml.D^D + A'^{H^H){Q^Q) + Al{H^Q){Q^H) 

+ {Al + hf)\e,,WQ^' + (/i?^ + gf,'')QlU:Q,pUp 

+ gf^''U^U:Q*^Q,, + gf/'D^D;Q*^Q,^ 

+ gl''D^DlQ]^Q,, + gl''U^U;D;D^ 

+ gl^'U.W^D^D, + Ai(Qtg)2 + \o,{D^Df + XqMQ^? 



+ \9lC^C^{U*y\^\^U\ (1) 



^ The effects of non-zero squark mixing in the case of a relatively light left squark doublet can complicate the 
calculation considerably. For a heavy left squark doublet at large values of the mixing parameters, the two-stage 
scenario is not realized |l5|. However, this effect may not persist for lower values of the mass of the left doublet as 
the contribution to the thermal mass of the scalars is changed, see HJ, |2^, On the other hand, as was noticed 
in previous studies of the phase transition, a non-zero value of the mixing parameters always weakens the strength 
of the transition. So lower values of the Higgs mass or the right handed stop mass are necessary to enhance the 
strength of the transition for non-zero squark mixing. Thus our results give upper bounds on the scalar masses. 
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Here H is the Higgs doublet field, U{D) is the right handed stop(sbottom) field, and Q is the third 
generation left squark doublet field. The longitudinal components of the SU(2) and SU(3) gauge 
fields are denoted by Ao and Co, respectively. The latin (greek) indices indicate SU(2) (SU(3)) 
components. As usual, the fields in eq. (|l]) are the static components of the scalar fields properly 
renormalized, the dimension of the fields in 3D is [GeV]^/^. Quartic couplings are of order gf{hf)T, 
having dimensions of [GeV]; here gi{ht) denotes a gauge ( top Yukawa) coupling. In the following 
we have not included the correction to the quartic coupling between the doublet Higgs field and 
the triplet scalar field Aq, or the corresponding correction for the SU(3) counterparts. We work 
throughout in the Landau gauge. For an analysis of the gauge dependence, we refer the reader to 
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We present the full relations between 3D coupling constants and masses in terms of the under- 
lying 4D parameters and the temperature. Partial results for the MSSM in the large-m^ limit can 
be found in ref. |]14| 0. 



2.1.1 Mass terms 

For the Higgs doublet we haveQ 



"^-3 = <(i + i^'T^-3^'T^) 

+ T^{^ + ^qOI + ^q9" + -h] + \{2hl sin^ /3 + 2A3 + A4)) 

~ Y^(6A< + 3K + m?,)/i^in2/3), (2) 

where the Higgs mass parameter is denoted by mn, and A = ^^-^^^^ ' cos 2/?, A3 = ^,A4 = — 
Lb = 21og^^ ~ 2 log y Q^gy , Lf = Lb + 4 log 2. Here JI is the mass scale defined by the modified 

^If a light higgsino is included, there will be important effects that are proportional to the top Yukawa coupling 
in the dimensionally reduced theory |2^. 

''We mostly neglect the hypercharge coupling g', throughout the paper. The only exception is in the contribution 
to the tree-level expression of the Higgs self-coupling A, as this latter quantity is fundamental in determining the 
strength of the phase transition. 
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minimal substraction scheme (MS) scheme. Similarly, for the third generation squark mass terms 
we have 



m 



Us 



U /4 



y 3 " 6 
1 « n ialml + 2hl sin^ /3(m^ + 



(3) 



— ql + — + -(?? + —hf(l + sin^ 3) 
16 2 9 ' 12 *^ ^ 



{^atm^Q + 6Aim^ + hfmfj + sin^ /^m^^ 



2^2 



(4) 



m^3 = my 1 + 4(7.? 



^6 



where the soft SUSY-breaking masses for the third generation left squark doublet and the right 

2 2 

handed sbottom are denoted by niq and m^, respectively, and \u = \i = The longitudinal 
components of the SU(2) and SU(3) gauge fields acquire thermal masses given by 
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(6) 



m: 



respectively, where Nsw = 4, Ns. 



Co - 

4,iV/ 



6 14 . 



ys \ 12 6 



(7) 



2.1.2 Couplings 

The 3D gauge coupling expressions for a dimensionally reduced SU (N) gauge theory can be found 



in ref. |]I4[. We include them for completeness 



9l 



Tg'iJi) 



487r2 



(22N - N,)^ ~ NfLf + N 



(8) 
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Here Ns is the number of scalar fields in the fundamental representation and Nf is the number of 
fermions. 

For the scalar quartic self-couplings, we have the following relations arising from the diagrams 
that have been shown in refs. ^ 0, 



Lb 



^ 2 Lb ru'i ■ 2 f3 

2 lovr^ 



167r2 



- T 



1 « 2 \ 77^91 + 12A^ + 3(A + A3A4 + Xj cos^ (3 + Xi sin* /3 



+ hi sin' /9(A3 + A4 sin' p) + ht sin* + - 



3 gt 



8 lovr^ 



167r2 



(9) 



A 



Us 



9 s rp 



8g! 



Lb 



167r2 

+ ht + ht sin' p) + ^gt 



T 



167r2 



23 4 13 4 2 2 4 

— <?« H 9« <7c 

36^' 12^" 3 



(10) 



73 



- T 



Lb 



-hi sin' (3gl + 2/1* sin* /3 + QXh] sin' 
3 



167r' 

+ /ii'(2A3 + A4 + /i?sin'/3 



(11) 



A 



Q 



X.T 1 



T 



:9 



2 Lb 

'167^2 



— 3/i' sin' 



Li 



167r^ 



IGtt^ 



+ 2A1A4 + A^ cos* + sin' /3 + 2A/i? sin' /3 + 2Ai/i? sin' /3 
+ /i* sin* P + 2hlXi sin* P + Xl sin* /^^ - -^2/i* sin* p] , 

/ lOTT^ J 



(12) 



^For contributions arising from the rest of the supersymmetric particles and the inclusion of Yukawa couplings 
for the other (s)quarks see ref. Eq]. 
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- T — ^ f 2AA4 cos^ p + 2X1X4 cos^ p + glXi cos^ /3 + 2A^ cos^ (3 

4 . \ 

+ AX^iXi cos^ (3 + -glXiCos^ [5 y (13) 



- T;TP^it:9Hhl + A4) sin^ /? + 2AA4 W /3 + 2A1A4 sin^ ^ + 2\l sin* /3 

/ lOTT^ \3 

+ 2/i^ sin^ 13 + 2A/i,2 sin^ /? + 2Xih^t sin^ /3 
+ h'^t sin^ /3 + 4/1^ A4 sin^ (3 + ^^(/i^ + A4) sin^ (3 + 2A^ sin^ p 

+ 4A3A4sin2/? + 4/ij2A3sm2/3), (14) 

hf ^ Ajs,V,T(l4,^^-3.?sin'^ + 4,I^) 

-T^TT^inOWt + ^4) sin^ + 2AA4 sin^ /3 + 2A1A4 sin^ + 2A^ sin^ /3 

2 lOTT^ Vo 

, 2/i^sin^/3 + 2A/i^sin2/? + 2Ai/it2sm^/3 
+ hi sin^ /3 + AhlXi sin^ /3 + ^^(/i^ + A4) sin^ /? + 2A^ sin^ /? 
+ 4 A3 A4 sin^ (3 + AhlX^ sin^ . 

We include the relations for the couplings among the heavy fields which will be integrated out 
at the second stage. These relations are needed only when the 2-loop contribution from these fields 
to the effective potential are included, see section 2.4. We obtain 



+ 



(15) 
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- + ^(,4 ^ 2hl\^ sin^ /? + /i2A4 sin^ (3 cos^ 2/? 

+ ^«.* + '';si"*f + ''?A4Sm*/3)+li^), (18) 

- l£T{l . 8»f Ji,) - T{^,{-hh! + + ^^). (19) 

- I,^t + l9l(2h^.-9l))+l4,). (21) 
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Ai 



+ 



23 



\9l + ^gl 



Lb \ 

(16^2) j 



2 V(16^2) 



hi 





u2 2 



j9t + Y^gt + -^glgl 



2\l + 2A3A4 + cos^ p + 2/1^2 A3 sin^ (5 



+ ht sin' P + 2hl\, sin' /3 + A^ sin^ (3 + ^gt) + ^gt 



(23) 



T 1 + 



+ 



:gl + Sfs 



91' 
6 



(167r2 



2 V(167r2 



2 
-/ 
3 



+ 2/i^A4 sin^ P+Xl sin' /? + ^^,^) + f^gi) 



2X1 + 2A3A4 + A^ cos^ + 2/1^ A3 sin^ /3 + /i^ sin^ 
13 A 13 



(24) 



fT(l+(?pi + 89,^)^ 

11 



V(16^2)^ t 



2„2 



Kg. 



+ TT^i's + 7f s + A4 cos^ (3 - 2/1; A4 sin^ /3 cos^ (3 - 2A| sin^ (3 cos' f3 + h; 



12 



4 



9: 



+ 2hl\i sin' (3 + \l sin^ (3 + ^.^^^ - j^gt^j + ^ (157^2); ' 



-Y 



(25) 



+ ^gl + 2A^ + 2A3A4 + 2/1^2 A3 sin^ /3 + 2hl\i sin^ (3 cos^ 

+ 2\lsin'Pcos'P+^-g'^ - \glgl) + ^^J^)- (26) 

A few technical comments are in order. We point out that the full one-loop contribution to the 
quartic coupling \eijH'^Q^\^ is given by the sum of eqs. (|1^) and (p^. Similarly, for the quartic 
coupling Q*i^U*Qi^Uj3 the full contributions arises from eqs. (|TB[) and (P^). It is also important 
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to note that there are off-diagonal (in colour space) gluonic contributions to the quartic couplings 
involving the strong gauge coupling. 

Additionally, when the full supersymmetric spectrum is not included then the running of each 
of the strong quartic couplings given in eqs. ([T6|)-(^6D is different. Although the gluino contribution 
is decoupled under our assumptions, we now write as a check the gluino contributions to the 
logarithmic part of the quartic couplings. This shows how, if one includes the full spectrum, then 
the relation between the beta-function coefficients is the same as the relation between the couplings. 

The gluino contribution, in units of T^^^g^, to the self-couplings is y. Using the same units, 
for gf^, gf^^ the contribution is —4, for gf^, gfj^ the contribution is y , for g^^^, g^^ the additional 
term is 14, and for g]!^ ,g'^^ the contribution is |. 

2.2 Second Stage 

Another simplification of the effective theory can be obtained by integrating out the scalar fields 
which are massive at the transition point. As we have seen the static modes corresponding to the 
scalar fields Q, D, Ao, Co, acquired thermal masses proportional to ~ gw{s)T, as a consequence of 
the integration out of the non-zero Matsubara modes. The second stage proceeds in exactly the 



same way as in reference |jT^. We include the additional corrections arising from the couplings we 
have considered. 

2.2.1 Couplings 

The final expression for the tree level 3D potential is given by 

VsD = m%H^H + \H,{H^Hf + mL^U^U + \u,{U^Uf + %H^HU^U, (27) 
where the scalar couplings are now 

16 oTrmAo oTcmQ.^ \ 

+ l{m'+iKr)+hi^^'i+h^K+lihi:r), (28) 
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W3 



Ac/3 ~ 



13 g. 



36 Svirmco 



UD\2 
S2 > 



{h?y - 2h?^gf^^ + 2h?^gi^ + + g^^ + 2{g, 



UD\2 

S2 I 



(29) 



(30) 



The 3D gauge couplings which appear in the SU(2) and SU(3) covariant derivatives of the 
effective theory are 



—2 _ 2/1 

9wz \ 



g 



9% = 9% ( 1 



9t 



gl 



g. 



S3 



IQnmco 247rmQ3 ASnniD^ 
2.2.2 Mass terms 

The one-loop contribution to the mass terms can be obtained directly as shown in ref. pi 



(31) 
(32) 



m 



Hi 



mrr 

167r 



(33) 



Us 



37r^ 



-(2<^ + 6^y^^)m^3. 



(34) 



Until now, our procedure has been exactly the same as in previous 3D reductions of the MSSM. 
However, in order to precisely fix the scales of the couplings that appear in the thermal polariza- 
tions of eqs. (2)- (3), one needs to perform a 2- loop evaluation of the effective potential. In their 
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calculation, Bodeker et al. took the values of the couplings in the screening parts of the 3D 



masses to be equal to the 3D values of these couplings. This is the correct result when two-loop 
corrections are included |2^, see eqs. ( p7| ) and ( |68D below. In addition, the mass parameters 
are renormalized in the 3D theory, 

m^3 (/i) = + JY^hm^ log ^ ' (35) 

(/i) = + .^^i-^/2„^ log (36) 

The expressions for the 2-loop beta functions f2mH > f2mu mass parameters have been given 

in ref. [^. As mentioned there in order to fix the values of the parameters Ah^ and A^/g we 



must employ the 2-loop effective potential of the 4D theory. In refs. [0, an estimate of 
~ A(73 ~ 7T was used. The expressions for the 2-loop effective potential in a iJ((/)-direction) 
and ?7(x-direction) background have been given for 4D in the paper by Carena et al. |T^, and for 



3D by Bodeker et al. In sections 2.3 and 2.4 we perform a two-loop calculation of the effective 
potential, incorporating all of the corrections to the 3D couplings obtained in the previous sections, 
to determine the exact values of A^-^ and A^g. 

We will analyse the effect on the critical temperatures when these corrections are included. 
Qualitatively we can say that if the net effect increases the value of Ah-^, then the critical temperature 
in the 0-direction decreases, and vice versa. A similar effect occurs in the ^-direction. Thus, as 
the range of values of m^^, which has been determined to give rise to a two-stage phase transition, 
is small, a more precise evaluation of the critical temperatures is of interest. This could have the 
effect of either reducing or enhancing the allowed range of values of the right stop mass, mi^, for 
which a two-stage phase transition can occur. 

2.3 Two-loop contributions 

The strategy we employ follows that of ref. . The idea is to use the 4D 2-loop effective potential 
in order to fix the scales in the 3D theory, and to use the 3D effective potential expressions for the 
Higgs and stop fields given in ref. [O] to analyse the phase transition. We calculate the unresummed 
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2-loop effective potential in order to include all 4D corrections to the mass parameters, resummation 
is automatically included in the calculation of the 2-loop effective potential in the 3D theory. We 
must also include the contributions to the 2-loop effective potential of the static modes, which have 
been integrated out at the second stage (includes the effects of resummation of the heavy fields). 

There are several effects that must be considered in order to obtain all of the contributions 
(constant and logarithmic) to the mass parameters. From the 4D effective potential one finds 
the two-loop contributions from the gauge bosons, Higgs, right handed stop, left handed squark 
doublet, right handed sbottom, and top quark. The expression for the 2-loop effective potential 



can be found in refs. \L4\ , |T^. In particular, within our approximations for the (^-direction. 



the appropriate expression is that given in ref. [0], because we are including the 3rd generation 



squark doublet and the right handed sbottom in the thermal bathQ. The main difference is that 
the D functions appearing below correspond to the unresummed expressions. Additionally we must 
include the effects arising at the second stage of reduction from the left handed squark doublet, the 
right handed sbottom, the scalar triplet and the scalar octet. We now derive the effective potential 
at finite temperature using the background fields and x = t^aU"', where we have chosen the unit 
vector in colour space tx" = (1, 0, 0). We now write the expressions in the shifted theory of the mass 
spectrum after the first stage of integration. The gauge bosons masses are 

2 ^ 2 i2 2 ^22 — 2 ^ 2 ^o^\ 

With no mixing in the Higgs sector, the Goldstone bosons and Higgs masses are 



2 2 , o \ a2 , u2 ■ 2 nX 

= mjj^ + 3A0 + sm p 



2 

2 2 2 I \ 2 I 7,2 • 2 n-r 

ml = ml + 3Xx' + hlsm-'P^. (St 



k2 



'For most of our analysis we will include the third-generation left squark doublet in the thermal bath. 
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The masses of the rest of the scalars contributing to the effective potential are given by 

ml^ =ml + {hi sin^ /3 + A3 + A4 sin^ + [h] - |) y , (39) 

<^ 3 =ml + {hi sin' /3 + A3 + A4 sin^ Z^) y + (f ) (40) 



m^^ = ,4 + (\, + M cos' 13)^+ [hi - I j ^, (41) 
2.3.1 (/)-direction 

The contributions to the terms in 'mjj,^{fi) proportional to g^, gl^Xn-j, and A^^ from Standard Model 
particles have been calculated in the paper by Kajantie et al. [^. The correct expression is 
obtained from eq. (151) of ||2T|, substituting A = \{g1, + g''^)cos2j3 and gy = /i^sin^/3. For the 



Standard Model, this already includes the finite contributions from counterterms. 

The additional corrections that arise from supersymmetric particles can be calculated using 
the 2-loop unresummed potential. The expressions in the integral form have been given in ref. 



1^ for zero squark mixing, and we include them for completeness. Our notation for the D- 



functions corresponds to that of ref. [^. The contributions from the 2-loop graphs containing 
supersymmetric particles are given below. For the 0-direction, we can drop the colour index of the 
squark masses: 
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(SSV) = — ^Nc[Dssv{'mt^,mL^'^w) + Dssv{'mi^,'mi^,mw) + 4:Dssv{mtL,mi^,mw)] 



- ^{N'^ - l)[Dssv{mi^,mi^,0) + Dssvim-,^,m-,^,0) 
+ Dssv{mi^,mi^,0) + Dssv{mi^, m^^, 0)], 



(45) 



(SSS) 



9^ \ 2 / ^2 ^ 2 



+ (/ij sin^ pfDsss{.min , '^Ir > i^h) + ( sin^ ^ + ^ cos 2/? ) Dsssi.miL ' "^6^. ' "^^^ ) 



-,(46) 



(5F) = --^2(^2_^)^^^^(^_^^o)_^^^^^^_^^Q)^^^^(^_^^Q)_^^^^(^_^ 



^glNj^Dsvi'miL^'^w) + ^5y("^6^,"^iy)], 



(47) 



(55) = ^iV,(2-iV,)Dc,5(m,-^,m-,J + /i?iV,[D55(m,-^,m,-J + D5s(mj^,m,-J] 

2 2 

+ (I" + j)Nc{N, + l)[Dss{mi^,raiJ + Dss{rn-,^.rn~J] 



6 



iV,(iV, + l)[Dss{mi^,miJ + Dss{m-,^,m-,J] 



+ Nc Q/ij sin^ (3 + ^gl cos 2/3^ [Dss{miL ' "^/i) + '^Dssi.mi^ , m^) + DssimL ' "^^)] 
- ^^cfi-^ cos2(3[Dss{mi^,mh) + 2Dss{mi^,m^) + Dss{mi^,m^)] 



+ -j^NcK sin (3[Dss{min^'mh) + 3L's5(mt-^, m^)]. 



(48) 



There are also counterterm contributions to the mass terms; for the 0-direction they correspond 
to eq. (B.3) of ref. P . 



SV = -^'^{?.gt + -g"). 
167r2 96^ 9 ^ ^ 



(49) 
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2.3.2 x-direction 

The two-loop unresummed effective potential in the x-direction is given by the following contri- 
butions^'!: 

(VVV) = -g^^^[iN,-2)Dvvv{mG,mG,0) + Dvvv{mG,mG,mG)], (50) 



ivvV) = -g',^[2{N, - l)D,,v{0,0,mG) + D,,y(0, 0, mc)], (51) 



(VV) = -gl^\2{N, - 2)Dyy(0,mG) + 2Dvv{mG,mG) + {N, - l)Dvy(mG, mc)], (52) 



c,2 

{SSV) = -YlDssvimtLi^'^hi^'^) + Dssv{mi^^,rnt^^,0) + Dssv{mi^,,mi^,,0) 
+ Dssvimi^, , mi^, , 0) + Dssv{mi^^ , m^^^ , 0) + Dssv{mi^., , mi^^ , 0) 

- 9s^[{Nc - l)Dssvirn^,rn^,rnG) + (N^ - l)Dssvim^,mu,mG) 

^ Dssv(rn^,mu,mG) + ■^Dssv{'m^,'m^,rnG) 
N,{N,-2)Dssv{m^,m^,0) + 2{N,-l)Dssv{mi^,,mi^,,mG) 

——Dssvirai^, , m,-^^ , m^) + ^^5sy ("^t^^ , , "^g) 
+ NlN,-2)Dssv{vfii^,.'mi^,.^) 

iVc - 1 _ 
+ 2{N^-\)DssY{mi^^,mi^^,mG) + Dssv{mi^^,mi^^,mG) 

+ ^^5sy(m^^2,m^^2,mG) + A^c(A^c - 2)L)s5v("^b^2'"^bL2''^) 



+ 
+ 
+ 



^ ' c 

^The expression given for the 4D effective potential would correspond to the usual resummed 2-loop 4D effective 
potential if we use the resummed expressions in the D functions appearing below. 

® As Tfii^^ = Tij^,^ , mi^^ = m-j^^^ , mi^^ — "i^^^ in the x-direction, we just multiply by a factor of 2 the contributions 
from these fields in some of the following expressions. 
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iVc - 1 

■ 2(7V, - l)Dssv{mi^,,mi^^,mG) + Dssv{mi^^,mi^^,mG) 



rrr (N — Ir 
(VVS) = -gl-^[{Nc-l)Dvvs{mG,mG,mu) + 2 — Dvvs{rnG,mG,mu) 

■ Nc{Nc-2)Dvvs{0,mG,m^) H ^ Dvvs{mG,mG,m^)], 



{SV) = -| 
1 



-^[2N,{N, - 2)Dsv{m^, 0) + {N, - l)[3Dsv{m^, mc) + Dsvim^, mc)] 

+ 2iV,(7V, - 2)Dsv{mi^„0) + {N, - l)[2Dsv{mi^,,mG) + 2Dsv{mi^„mG)] 
+ ['^^svimi^, , rfiG) + 2Dsv{mi^, , mo)] 

+ 2N,{N, - 2)Dsv{mi^,,0) + (AT, - l)[2L>5y(m^^^, mc) + 2L>5y(m^^^, mc)] 

+ 2Nc{Nc - 2)Dsv{mi^,,0) + {Nc - l)[2Dsv{mi^,,mG) + 2Dsv{mi^,,mG)] 
+ , mo) + 2Dsvimi^, , mo)] 

+ Dsvimi^^ , 0) + Dsvimi^, , 0) + Dsv{mi^, , 0)] , 



(SSS) = -A^x i^Dsssimu, ruu, rriu) + (2iVc - l)Dsss{mu, m^, m^)] 
- X^i-^h^ sin^ p[Dsss{mu,mh,mh) + 3Dsss{mu,m^,m^)] 
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+ 2(/it - -glfDsss{m^, mi^^,mi^J + 2{-^fDsss{mu, m^^^, m^^J 

1 (,2 

+ 2{^fDsss{m^:m-,^^,m-,Ji (56) 

^SS) = ll^^ + ^y8Dss{mi^,,miJ + 24Dss{mi^,,miJ + 16Dss{mi^,,m 
+ SDss{mi^,,miJ + 2ADss{mi^,,mt,J + 16D55(mg^^, mg^J] 

+ \i^)[8Dssimi^,,miJ + 2ADssimi^,,miJ + IGDssim^^.^miJ 

+ h?t [Dss{mi^^ , mi^^ ) + Dss{mi:^^ , ) + Dss {mi^^ , mi^ ) ] 
+ h\ {Dss{mi^^ , mi^^ ) + Dss{mi^^ , mi^^ ) + ^55 ("^^^3 , ^,-^3 ) ] 

+ (2^* + 8^"' cos 2/?^ [i:>S5(m4^^, m^) + Dss{mi^^,mh) + Dss{mi^^,mh) 

+ 2{Dss{mi^-^,m^) + Dssim-^^^,m^) + Dssimi^^,m^)) 
) + Dssimi^^.m^) + i:'55(m£^3, m^)] 

- Qi/^ cos 2/3^ [i:)55(m5^^,m^) + Dss{mi^^,mh) + Dss{mi^^,mh) 
+ 2{Dss{mi^^ , m^) + Dss{mtL2 ' "^'t) + Dssimf^^ , m^)) 

+ ^55 ( 

) + Dss{mtR2^^^) + Dssimins^^^))] + ~ ^c)[-D55("7t^i, "^fe^J 
+ Dss{mi^^,mi^^) + Dss{mi^^,miJ]. (57) 
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The counterterm contribution to the mass term comes from an analogous contribution to that of 
eq. (|9D and isf\ 

y2 146 , 

2.4 Integration over the heavy scale 



The second part of the calculation arises, as noticed in the paper by Kajantie et al. ||2T[|: when the 
"heavy" particles have been integrated out their contributions to the 3D mass parameters should 
also be included, as they can substantially vary the value of the parameters Ah^jAu^. In order to do 
this we must calculate the 2-loop contributions to the effective potential in the 0- and ^-directions 
from the heavy fields: Q, D, Co, Ao 

For the (;/)-direction the expression of the effective potential to 2-loops from the heavy particles 
can be deduced from the expressions in the paper by Espinosa . For the x-direction this is new. 
The masses in the shifted theory are now given by 



<, = -ga + (^f + A3 + A:)f + {hT + 9?!" + i?r)y ^ (59) 



= + (hf + As + ADy + (^7?,'^) y, (60) 
= + (A3 + Al)^ + [h^ + 9?!" + ^r)y , (61) 
<.3 = + (As + ADy + (^72^)y , (62) 



rrir = m 



(<" + C")^' (63) 



bill ~ -03 ^ vysi ^ J 2 



<.,,=^% + (92r)\- (64) 



There is an additional counterterm contribution in the ^-direction for the case of a Hght higgsino. 
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The expressions for the rest of the fields are given in . The 2-loop contributions from the heavy 
scale are given below. We stress that the D-integrals in eqs. ( |65D and ( p6| ) are just 3D integrals, 
our notation follows that of refs. 



2.4.1 0-direction 

3 

yheavy ^ -^^ j^^^^^^^^ _^ SOLsim^, "l^o)] + '^gl-.D Lvimw , ^Aq) 

3 3 

~ j^9ls(p'^DLLs{'mh,mAo,mAo) - -gl^DLLv{mw,mA„,mAo) 

gl 

- -^Nc[Dssv{mtL^mL^^w) + Dssv{mi^,mi^,mw) + 4:Dssv{miL^^bL^^w)] 

- -f{Nl - l){Dssv{mi^.mi^,Q) + Dssv{mi^,m-,^,Q) + Dssv{m~,^. m-,^, 0)] 

- [(/if + As + AlfDsssirrii^ , rrii^ , mn) + (A3 + AiyDsssimj^^ , mj^^ , m^) 
+ {h^ - A2 + A^)2Z}55s(mi^,m^^,mO]^Ar, 

3 

- gfi'^3^c[I^5v("^t^, "iw) + Dsvi'mi^^'mw)] 

+ (2Ai + + gl'i)N,{2 - N,)Dssimi^,m-J 

+ {h?'' + gf!" + 3gf^'')N,[Dssimi^ , m,-J + Dssim-,^ , m^J] 

+ (Ai + AQjA^,(A^e + l)[Dss{mi^,miJ + D5s(m^^,m^J] 

+ XD:MN, + l)[Dss{mi^,m-J] 

+ Ar,i(/if' + A3 + A^)[D55(m,-^, m;,) + 2Dss{mi^,m^) + Ds5(mj^, m^)] 

+ ^N^{A3 + Al)[Dss{mi^,mh) + 2Dss{mi^,m^) + Dss{mi^,m^)] 
+ {g?,'' + 3g?J')N,[Dss{mi^ ,m~J + Dss{m~,, , m~J] 

+ {gl"" + 3gl'')K[Dss{mi^.m-,^)]. (65) 
^°Our convention for the functions Dyvv, Dyvs is that of fl^ . 
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2.4.2 ^'direction 



■^heavy 



+ 



+ 



+ 



-j-UNc- l)[DLs{mu,rnco) + DLs{rnu,,'mco)] 
■^[4DLs{m^,mco) + 2DLs{mu,mco)] 



- 9: 



2 Nc 



(AT, -1)2 

2 I (-^c - i)DLLs{'mco,mco,mu) + 2 — DLLs{mco,mco,my,) 

(Nc - 2)2 

+ Nc{Nc - 2)Dlls{0, mca.m^) H "-j^ DLLs{rnco,mca,mS) 

- 9s3-Y^-DLv{mco,rnG) - DLv{mco,mG) - {Nc - l)DLv{mco,mG)] 

- 9%'^[DLLv{'mco,mc,,,mG) + 2DLLvimco,mco, ruG) 
+ DLLvirnCf,,mco,0) + 2DLLv{'mco,0,mG)] 

9^ 

- -f- [Dssv {rrit^, , rrif^^ , 0) + Dssv{mt^^ , mi^^ , 0) + Dssv{mi^^ , mi^^ , 0) 
+ Dssv{mi^^,rni^,,0) + Dssv{mi^^, m^,^^, 0) + Dssv{mi^,, m^,^,, 0) 



+ Dssv {rrii^, , ,rnG) + Nc{Nc - 2) Dssv {rrit^, , "if^^ , 0) 



+2{Nc - l)Dssv{mf^^ , mi^^.niG) + 



+ 2{Nc - l)Dssv{mi^^,mi^^, mc) + ^'^^ ^ Dssv{mi^, , mi^, , mo) 

0) + 2{Nc - l)Dssv{mi^,,mi^^,mG) 



+ —Dssv{m~,^^,m-,^^,mG) 
+ Nc{Nc-2)Dssv{mi^,,mi^^ 



N. 



Dssv{mi^^,mi^^,mG) 



Jj-Dssv{mi^,,mi^,,mG) + Nc{Nc - 2)Dssv{mi^,,mi^,, 0) 



- Y^ih?"" + 9?!" + 9%''fDsss{m^, ^i^.^^ij + 2(/^r + ^?,^)^i5555(m., m,^^ , m.^J 
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fl2 r 

- ^^2iVe(iVc-2)D5yK^,,0) + (iVe- l)[2D5yK-^,,mG) +2I)5KK^„mG)] 
+ ^[2-D5y(rnf^2,^G) + 2Dsy(mt^^,mG)] 

+ 2iV,(iV, - 2)D5v(m^^^,0) + (iV, - \){2Dsv{m-^^^,mG) + 2I)5y(m^^^, mc)] 

+ 2Ar,(iVe - 1)Dsv{mi^^,^) + (TVe - l)[2L»5y(mg^^, mc) + 2D5y(m^^^, mc)] 

- ^^i[^5y(mf,^, 0) + L'5y(mi,,, 0) + ^5^(771,-^3, 0) 
+ i^sy (777^^^, 0) + D5vK^,,0) + L'5y (777^^3,0)] 

+ ^(Ai + \Q,){^Dss{mi^,.vaiJ + 2^Dss{mi^,.raiJ + 16^55(777,-^,, 777,-^ J 

+ 8^55(777^^^, 777g^J + 24i:'5s(r77g^^, 777^^J + 16^55(777^^^ , 777-^^ J] 

+ -^{^D^)^^Ss{ml^^,ml^^ + 24D55(777g^^,r77g^J + 16^55(7775^^, 7775^J] 

+ {^T + gfWssirni^^.miJ + ^55(777,-^^, Tn.-^J + ^55(^,^3, 7^,-^3)] 

+ [^5S (^t^i , mi^^ ) + i^55 (m,-^! , mi^^ ) + ^55 {mi^^ , 777,-^3 ) 

+ ^55 (^77,-^^ , 777,-^^ ) + D ss (mf^, , mf^^ ) + Dss (mf^, , vai^^ ) 

+ ^55 (r77t^3 , 777,-^^ ) + DsS {vdl^^ , 777,'^^ ) + D55 (777,-^3 , 777,^3 ) 

+ ^ss {mi^^ , 777i^^ ) + Dss {mi^, , m,-^^ ) + Dss {mi^, , ^,-^3 ) 

+ Dss {rrii^, , mf^^ ) + Dss (m^^^ , 777,-^^ ) + Dss {mi^, , m^^ ) + Dss {mi^, , rui^^ ) 

+ Dss{mi^,,miJ + ^^55(^5^3, 777,-^3)] 
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+ 2(Ds5( 

) + -Ds5(W-ti3,"^7r)] 

+ Dssi 

+ 2Ai(2 - N,)[Dssimi^„m-,J + Z^^sK^^' ^^fe.^) + Dssimi^,,m-,J] 

+ Dss ' "^fefli ) + Dss {rrii^, , mj^^, ) + ^55 (m^^^ , m^^^ ) 
+ {rriir., , "^fe^, ) + Dss (rrii^, , "^6^, ) 

+ Dss {mi^, , m^^^ ) + Dss (mj^^, , mj^^, ) + ^55 ('^b^, , m-^^, ) + ("^6^, , mj^^, ) 

+ I^5s(m^,,, m^^,) + Dssim-,^^,m-,J + Dss{m-,^^,miJ + Dss{mi^^,m-,J 

+ dl," [Dss {mini ' ^bn. ) + ("^?fl2 ' "^6«2 ) + (^^^3 , mi J] 
+ ^s^,^ [Dss {mi^, , m^^^ ) + Dss {mi^, , "^6^, ) + Dss {mm, ' "^l^, ) 
+ ' "^bfli ) + {m-t^^ , mj,^^ ) + Dss {mi^^ , mj,^^ ) 

+ Dss {m-t^^ , mi^^ ) + Dss {mi^.^ , m^^^ ) 

+ Dss{mi^^^,mi^J\. (66) 
2.4.3 Mass terms 

Using the results presented in the previous sections, we can finally write the full expressions for the 
mass terms of eqs. (|35D and (^) 
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T2 / ,/146 2 11 



(167r2) V ' V 27 3 54 

1 / 1 4 /29 / / 3T \ 1\ A / 3T \ 1 

+ (i6^^(-4*Hyr<(2;;^j) +'^+ 2) +<'°<s^r) 2 



1 

+ c+ - 



1 



+ ^(^4-4e)-24e)> (67) 



+ l/i/ + l(/if + 2Af + A^ + A^ + 73) 



+ (§^»3 + 9 Ah3^^3 - 12Ai3 - 37^ + 8gl%) (log 

/ 4 137 9, ^ 1, 3,2 



+ c 
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- 2^?,23(-4(/.f^ + + Al) - 4(A« + AD) (^log(^^^) + c + ^) 

- 3((/.,^ + A? + Al)2 + (A^ + Air + (hf -K + Kf) + 

+ (t£) (^^'^' ^ + + ^^t - ^^t) , (68) 

where 



hi 



htsin' (3T[ 1 



8 (16^2 



(l2hlsiii^P-Qgl 



64 



9l]Lf 



64 



+ c?^ (2 + 28 log 2) - 96A log 2 + 16/i; sin^ /? log 2 ^,^(4 log 2-3) 



9 



. 9 C(2) 
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-0.348725. 



(69) 
(70) 



The exact values of the parameters and Af/g in eqs. ( pSj) and (^), depend on the particle 
content of the theory and on the input parameters. In fact, the dependence on tan/? and mt^, 
for the range of values we are interested in, is weak. The dependence on mg and m/) is stronger. 
For niQ = = 300 GeV, the corresponding values of the parameters are in the range = 
(1.6 - 1.8)r, and Ku^ = (6.6 - 6.9)T. 

2.5 Zero- Temperature Renormalization 

In order to complete the matching of the 3D parameters to the 4D physical parameters, we must 
renormalize the zero-temperature theory. We will not go into the details of the renormalization, 
but refer the reader to the literature in which the pole masses for the relevant particles of our 
calculation have been obtained considering the full particle spectrum of the MSSM [0, |3^. We 
use the expressions given in ref. ||31|, keeping only the top Yukawa coupling, in the appropriate 
(large-m^) limit. We have kept all T^-dependent contributions to order gf{h^) but have neglected 
the constant contributions, which are not multiplied by hf or gg. In this way all explicit dependence 
on Jl is cancelled at one-loop when we relate the 3D parameters to pole masses. 
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Another zero-temperature constraint that we impose is the stabihty of the physical vacuum. 
In principle, a metastable region exists in which the colour-breaking minimum is lower than the 
physical one at zero temperature. If the time for the transition to this lower minimum is greater 
than the age of the Universe then this region of parameter space is also acceptable. However, we 
will not consider these issues in the present paper. The constraint for absolute stability can be 



obtained by studying the effective potential at zero temperature [[I^, This gives the constraint 
—rriij < {mfj)'^, where 

„= = f!!%M)"'. (71) 
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3 Results 



With the previous results we can now analyse the phase transition. In fig. Q we show the critical 
temperatures for the transitions in the 0- and x-directions as a function of the right handed stop 
pole mass fni^, for tan/? = 3,5,12. We find that, for mg ~ 300 GeV, there still is a region in 
which a two-stage phase transition can occur. This region is to the left of the crossing points of 
the curves. With respect to the work of ref. [I^ the structure of the phase diagram is preserved. 



although it is slightly shifted towards higher values of the right stop mass. Note that there is a 



considerable difference between our values of the critical temperatures and those in ref. []14|]. Our 
analysis concludes that the structure of the phase diagram is robust to small additional corrections. 
This structure is maintained also for itlq = mD = 1 TeVQ The total effect does not substantially 
increase or decrease the range of values of the right handed stop mass for which a two-stage phase 
transition can occur. However, the exact location of this small range in the value of mi^ depends 
on the value of the third-generation left handed squark doublet mass. In fig. |^ we give the values 
of |; for three different values of tan (3. As expected, the strength of the phase transition has a 
weak dependence on the values of the scales that have been fixed in our calculation, and only slight 



^^In this case the third-generation left handed squark doublet and right handed sbottom are decoupled from the 
thermal bath. The relations obtained from the dimensional-reduction procedure can be deduced from all of the 
formulae presented in the previous sections. All contributions from these fields are suppressed at finite temperature; 
however, a residual dependence on mq as a consequence of a zero-temperature effect persists for the scalar Higgs 
self couphng A 
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differences are observed with respect to previous analyses. 

Figure ^ shows hues of |; = 1 for three different cases in the rrih-mi^ plane. The phase transition 
is sufficiently strong for electroweak baryogenesis to the left of the solid (dotted) line for mq = 300 
GeV (mg = 1 TeV), using the results obtained in this paper. The dashed line is the result using 
the approximations of ref. [|l4l, for mq = 300 GeV. We can see that the full effect of the corrections 
we have included on the strength of the phase transition is small. The end-points of the lines 
correspond to the maximum value of the Higgs mass that is reached by the effect of the zero- 
temperature radiative corrections for a given value of niq, and the y-axis starts at the value of the 
experimental limit on the Higgs mass^ p3| . 

The allowed region in parameter space is shown in figs. 4 and 5, given the current experimental 
limits on the Higgs mass, for two different values of mq. The region on the left of the solid line 
indicates when a sufficiently strong first-order phase transition occurs. The dotted line gives the 
condition for absolute stability of the physical vacuum. As explained above, to the left of this line 
the colour-breaking minimum is lower than the physical one at zero-temperature. The dashed line 
is obtained when the critical temperatures of the transitions in the (j)- and x-directions are the same. 
A two-stage phase transition occurs to the left of the dashed line. Note that, unlike the results of 
ref. |T^ for mq = 1 TeV and zero squark mixing, the dashed and solid lines do not intersect. 

As mentioned in the introduction, for a sufficiently heavy right handed stop field, an effective 
theory with a single light scalar Higgs doublet field can be constructed. In the appendix we give 
the formulae that modify the equations presented in the previous sections. There are two ways of 
constructing this effective theory. We can either integrate out the right handed stop field simultane- 
ously with the other heavy fields, or perform a third stage of reduction and integrate out the right 
stop field separately. The basic change in the equations of the appendix will be to replace the 3D 
coupling and masses by the 3D barred couplings and masses of section 2.2. It is clear, as mentioned 
above, that the results obtained with these effective potentials are unreliable for low values of the 
right handed stop mass, as perturbation theory is no longer under control. In terms of the 3D 
parameters we note in particular that the difference between Xh-^ and Xh^ is considerable, owing to 
the difference in mfj^ and mfj^. However, comparing the results obtained with these three separate 
-•^^As we are working in the large-mA limit, the Standard Model bound on the Higgs mass is used. 
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approximations, we can obtain a more precise value of the right handed stop mass, for which the 
effective theory with a single light scalar field is a valid description. Recall that the one-loop esti- 
mate gave as a lower limit the value ^177 GeV [0, In fig. 6 we present the plot of the 
ratio of the vacuum expectation value of the Higgs doublet to the temperature as a function of the 
right handed stop mass obtained for the three approximations for the effective potential mentioned 
above. In this plot we have taken tan (3 = 5. The solid line is the result obtained when using the 
2-loop effective potential derived in the main part of the paper. The dotted line corresponds to the 
results obtained after integrating out the right stop field simultaneously with the other heavy fields. 
The dashed line is the result obtained when a separate third stage of reduction is performed. We see 
that only for very large values of the right handed stop mass are the results basically the same. As 
we move towards lower values of m^^, the lines very quickly start to diverge. In particular, we see 
that for m^^ = 177 GeV {mu ~ 50 GeV) the results using the effective potentials for a single light 
scalar are completely unreliable. To understand why the results of the dashed and dotted lines are 
so different for smaller values of ttli^, we recall that the strength of the transition is dominated by 
the value of Xh^ in the final 3D effective theory for a single light scalar field at the phase transition. 
The value of is much smaller when a third stage of reduction is performed, as the value of m^^ 
is much smaller too. That is, the net effect artificially strengthens the first-order phase transition. 

4 Conclusions 

We have performed a full two-loop dimensional reduction of 4D MSSM parameters to the 3D 
couplings and masses of the effective theory. In this way, we have fixed the scales appearing in the 
3D mass terms that are due to the thermal polarizations and the super-renormalizability of the 3D 
theory. The values of the parameters A^^ and A^/g can vary significantly for different values of the 
input parameters and the particle content of the theory, thus modifying the critical temperatures 
of the transitions. We have compared our results with previous analyses. We conclude that the 
corrections relevant to the preservation of the baryon asymmetry are small and that the allowed 
range of masses is rrih ^ 110 GeV and m^^ < rrit, in complete agreement with previous results. We 
find that the phase structure diagram still allows a possible two-stage phase transition for a small 
range of values of . This range of values is shifted compared to previous results. However, 
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whether or not the transition actually occurs must be explicitly checked. Initial lattice analysis 
suggests that the second stage of the transition is extremely strong and thus this transition might 
not have taken place on cosmological time scales. Consequently, this region of parameter space for 
electroweak baryogenesis would be excluded. From the comparison of the results obtained using 
the different approximations to the effective potential we can conclude that for m^^ < 210 GeV, 
the results obtained for the strength of the phase transition with a single light Higgs field at the 
transition point are unreliable. 
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Appendix 

Integrating out the right handed stop 

The additional corrections that must be included in the case of integrating out the right handed 
stop are given below. The Higgs self-coupling is modified by 

A3 = As - t|-73'^. (72) 
The one-loop correction to the mass term is 

^3 = "^3 - ^73?^c/3- (73) 



The relevant contributions from the 2-loop graphs have been given in ref. |T2[ with the following 



modifications, which include the effects of the dimensional-reduction procedure. For the (SSV) and 
(SV) terms substitute Qs gs^, for the (SSS) term substitute /if sin^/3 — > 7I for the (SS) term 
substitute ^ Xu^ and hf sin^ (3 ^ j^. 



31 



References 

[1] V.A. Kuzmin, V.A. Rubakov, and M.E. Shaposhnikov. Phys. Lett, B155:36, 1985. 
[2] M.E. Shaposhnikov. Nucl. Phys., B287:757, 1987. 

[3] A. Cohen, D.B. Kaplan, and A.E. Nelson. Annu. Rev. Nucl. Part. Sci., 43:27, 1993. 
[4] V.A. Rubakov and M. Shaposhnikov. Usp. Fiz. Nauk, 166:493, 1996. 
[5] M. Quiros. Helv. Phys. Acta, 67:451, 1994. 
[6] S. Myint. Phys. Lett, B287:325, 1992. 
[7] G.F. Giudice. Phys. Rev., D45, 1992. 

[8] J.R. Espinosa, M. Quiros, and F. Zwirner. Phys. Lett, B307:106, 1993. 

[9] A. Brignole, J.R. Espinosa, M. Quiros, and F. Zwirner. Phys. Lett, B324:181, 1994. 
[10] M. Carena, M. Quiros, and C.E.M. Wagner. Phys. Lett, B380:81, 1996. 
[11] D. Delepine, J.-M. Gerard, R. Gonzalez Felipe, and J. Weyers. Phys. Lett, B386:183, 1996. 
[12] J.R. Espinosa. Nucl. Phys., B475:273, 1996. 
[13] B. de Carlos and J.R. Espinosa. Nucl. Phys., B503:24, 1997. 
[14] D. Bodeker, R John, M. Laine, and M.G. Schmidt. Nucl. Phys., B497:387, 1997. 
[15] M. Carena, M. Quiros, and C.E.M. Wagner, preprint ( |hep-ph/971040lD , 1997. 



[16] J. Cline and G.D. Moore, preprint (|hep-ph/980635^ ), 1998. 

[17] R Ginsparg. Nucl. Phys., B170:388, 1980. 

[18] T. Applequist and R. Pisarski. Phys. Rev., D23:2305, 1981. 



32 



[19] S. Nadkarni. Phys. Rev., D27:917, 1983. 

[20] K. Farakos, K. Kajantie, K. Rummukainen, and M.E. Shaposhnikov. Nucl. Phys., B425:67, 
1994. 

[21] K. Kajantie, M. Laine, K. Rummukainen, and M.E. Shaposhnikov. Nucl. Phys., B458:90, 1996. 

[22] A. Jakovac and A. Patkos. Phys. Lett, B334:334, 1994. 

[23] E. Braaten and A. Nieto. Phys. Rev., D51:6990, 1995. 

[24] M. Laine. Nucl. Phys., B481:43, 1996. 

[25] J. Chne and K. Kainulainen. Nucl. Phys., B482:73, 1996. 

[26] M. Losada. Phys. Rev., 56:2893, 1997. 

[27] G. Farrar and M. Losada. Phys. Lett, B406:60, 1997. 

[28] M. Losada. in preparation. 

[29] P. Arnold and O. Espinosa. Phys. Rev., D47:3546, 1993. 

[30] M. Laine and K. Rummukainen. preprint (|hep-ph/9804255| ), (|hep-lat/9804019| ), 1998. 



[31] D. Pierce, J. Bagger, K. Matchev, and R. Zhang. Nucl. Phys., B491:3, 1997. 
[32] A. Donini. Nucl. Phys., B467:3, 1996. 

[33] V. Ruhlmann-Kleider, DELPHI Collaboration. Contribution to the Moriond Conference on 
Electroweak Interactions, 1998. 

[ittp: / / delphi www.cern.ch/ delfigs/searches/ vanina980325/ mor98.ps 



33 



93.0 - 



88.0 - 



83.0 - 



78.0 - 




73.0 I — • — • — • — • — ' — • — • — • — • — ' — • — • — • — • — ^ 
157.0 162.0 167.0 172.0 



Figure 1: Critical temperatures in the 0- (solid) and (dotted) directions as functions of m^^ for 
tan (3 = 3,5, 12 and tuq = 300 GeV. 
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Figure 2: Plot of |; as a function of m^^ in the 0- (solid line) and (dotted line) directions for 
tan /? = 3, 5, 12 and ttiq = 300 GeV. For a given value of tan/3 the lines end at the same value of 
the right handed stop mass. 
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Figure 3: Contours of |; = 1 in the rrih-mi^ plane. The sohd (dotted) hne corresponds to the results 
obtained within our approximations for mq = 300 GeV (ITeV). The dashed line is the result using 
the approximations of ref. |T^ for mq = 300 GeV. The region to the left of the lines gives a 



sufficiently strong first-order phase transition, for a given value of mq. 
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Figure 4: Allowed region in nih-mi^ plane for mq = 300 GeV. To the left of the solid line there is a 
sufficiently strong first-order phase transition, to the right of the dotted line the physical vacuum is 
absolutely stable. The dashed line separates the region for which a two-stage phase transition can 
occur. 
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Figure 6: Plot of |; as a function of m^^ in the 0-direction for three cases. The sohd hne corresponds 
to the result obtained with the effective potential given in the main part of the paper. The dotted 
line corresponds to the results obtained after integrating out the right stop field simultaneously 
with the other heavy fields. The dashed line is the result obtained when a separate third stage of 
reduction is performed. 
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